Electronic structure and optical properties of lightweight metal hydrides 
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We study the electronic structures and dielectric functions of the simple hydrides LiH, NaH, MgH 2 
and AIH3, and the complex hydrides LiaAIHe, NaaAlHg, L1AIH4, NaAlHU and Mg(AlH 4 ) 2 , using 
first principles density functional theory and GW calculations. All these compounds are large gap 
insulators with GW single particle band gaps varying from 3.5 eV in AIH3 to 6.5 eV in the MAIH4 
compounds. The valence bands are dominated by the hydrogen atoms, whereas the conduction bands 
have mixed contributions from the hydrogens and the metal cations. The electronic structure of 
the aluminium compounds is determined mainly by aluminium hydride complexes and their mutual 
interactions. Despite considerable differences between the band structures and the band gaps of 
the various compounds, their optical responses are qualitatively similar. In most of the spectra the 
optical absorption rises sharply above 6 eV and has a strong peak around 8 eV. The quantitative 
differences in the optical spectra are interpreted in terms of the structure and the electronic structure 
of the compounds. 
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I. INTRODUCTION 

The large scale utilization of hydrogen as a fuel cru- 
cially depends on the development of compact hydrogen 
storage materials with a high mass content of hydrogen. 1 
Hydrides of group I-III metals in the upper rows of the 
periodic table could meet this requirement. These met- 
als are sufficiently light for their hydrides to have a large 
gravimetric hydrogen density; for instance, MgH 2 con- 
tains 7.7 wt.% hydrogen. One must be able to extract hy- 
drogen at a moderate temperature, however, and there- 
fore the metal hydride should be neither too stable nor 
too unstable. Simple metal hydrides do not satisfy this 
demand. For example, the binding energy of MgH 2 is too 
largej^ whereas the binding energy of AIH3 is close to 
zero^ 

This has stimulated research into binary intermetallic 
hydrides such as the alanates MAIH4, M'(A1H 4 ) 2 , with 
M and M' a light alkali and alkaline earth metal, respec- 
tively. Some of the properties of these compounds have 
indeed improved as compared to the simple hydrides, but 
the compound that meets both the stability and the stor- 
age capacity demands has not been found yet. Whereas 
sodium alanate, NaAlH.4, releases hydrogen in two reac- 
tion stages with enthalpies close to the ideal value, its 
active gravimetric hydrogen density is only 5.5 wt.%i£*& 
Magnesium alanate, Mg(AlH 4 ) 2 , and lithium alanate, 
LiAlH4, have a higher active gravimetric hydrogen den- 
sity of 7.0 and 8.0 wt.%, respectively. However, they are 
not sufficiently stable with respect to decomposition into 
simpler hydrides .LSuLifliiiiiSi 

A suitable ternary intermetallic hydride might satisfy 
all requirements. The number of possible ternary com- 
pounds is very large, however, and searching for the opti- 
mal composition becomes very tedious, unless one uses a 



combinatorial approach. Such a technique has been pro- 
posed recently, in which thin films are grown with tun- 
able composition gradients^ The composition is then a 
function of the position on the film. This avoids having 
to synthesize all compositions individually, but one still 
needs to be able to identify the most promising ones. It is 
proposed that identification can be based upon the opti- 
cal properties of suitable metal hydrides being very differ- 
ent from those of their host metals^ This has first been 
demonstrated conclusively for YH^^l and since then for 
a number of other so-called "switchable mirror" rare 
earth and transition metal compounds ^iLiLi^iiSiSSiSi 
The compounds that absorb the maximum amount of 
hydrogen, become semiconductors or insulators. 

If one applies this technique to group I-III metal hy- 
drides, it is relevant to know how the optical properties of 
these compounds depend on their composition and struc- 
ture. In this paper we report a systematic first prin- 
ciples study of the band gaps, the electronic structure 
and the optical properties of group I-III metal hydrides. 
Band gaps and single particle excitations are calculated 
within the GW quasi-particle approach; optical excita- 
tions are obtained using the random phase approxima- 
tion (RPA). We focus upon a number of elements that 
are of interest for lightweight metal hydrides, i.e., Li, Na, 
Mg and Al. In particular, we consider the series of simple 
metal hydrides LiH, NaH, MgH 2 , AIH3 and the binary 
metal hydrides Li 3 AlH 6 , Na 3 AlH 6 , LiAlH 4 , NaAlH 4 and 
Mg(AlH 4 ) 2 , 6 i 8 i 22 i 23 i 24 i 25 i 26 i 27 i 28 The trends in the optical 
spectra and electronic structure are discussed and inter- 
preted in terms of the structure and bonding of the ma- 
terials. 

This paper is organized as follows. In Sec. II we out- 
line the computational methods used in our study. The 
results are presented in Sec. Ill, first for the simple hy- 
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drides, then for the binary M3AIH6 hydrides and finally 
for the MAIH4 alanates. Sees. IV and V contain the dis- 
cussion and a summary. 



II. COMPUTATIONAL METHODS 

The results discussed in this paper are obtained using 
a combination of density functional theory (DFT) and 
GW calculations. DFT is used at GGA level to optimize 
the ground state structure and obtain single particle wave 
functions to be used in the calculation of the optical re- 
sponse. GW is used to generate single particle excitation 
energies within the quasi-particle (QP) approximation, 
starting from DFT/LDA wavefunctions and eigenvalues. 
The optical response is given by the frequency dependent 
dielectric function, which is calculated within the random 
phase approximation (RPA). In the latter we use single 
particle wave functions and excitation energies and ne- 
glect exciton and local field effects. The main difference 
between the DFT and the QP excitation spectra is the 
size of the fundamental band gap between occupied and 
unoccupied states, whereas the dispersion of the bands 
is quite similar. We use a scissors operator to the DFT 
eigenvalues to approximate QP excitation energies on a 
dense grid in the Brillouin zone, which is required to cal- 
culate the dielectric function. 



A. DFT calculations 

DFT total energies are calculated with the PW91 gen- 
eralized gradient approximation (GGA) functional 2 ^ and 
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as 



the projector augmented wave (PAW) method 
implemented in the Vienna Ab initio Simulation Pack- 
age (VASP) i? 2 '??i? 4 We use standard frozen core PAW 
potentials and a plane wave basis set with a kinetic en- 
ergy cutoff of 312 eV. The tetrahedron scheme is applied 
for the Brillouin zone integration using k-point meshes 
with a spacing between 0.01 and 0.03 A^ 1 . The cell pa- 
rameters and the atomic positions within a unit cell are 
optimized by minimizing the total energy except for the 
alkali alanates, where only the atomic positions are opti- 
mized and the cell parameters are taken from experiment. 
These optimized structures are used as input for the GW 
calculations. 



B. GW calculations 

DFT calculations generally give good results for 
ground state properties, but not for excited states. The 
electronic band gap, for instance, can be underestimated 
by ~ 50%, and even more than that for small band gap 
materials. This stems from an unjustified interpretation 
of the DFT (Kohn-Sham) eigenvalues as single particle 
excitation energies. The latter can be obtained from 



the quasi-particle (QP) equation, which involves the non- 
local, energy dependent self-energy £, 



1. 
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-0nk(r) + 

y dr'E(r, r';enk)V'nk(r') = e„k^ n k(r), (1) 



where w cx t stands for the sum of all nuclear or ionic poten- 
tials, Vh is the electrostatic or Hartree potential resulting 
from the electrons, ?/>nk(r) is the QP wave function, and 
e„k is the QP energy, i.e., the single particle excitation 
energy. In practice Eq. is solved using a number of 
approximations. The GW technique approximates the 
self-energy £ by a dynamically screened exchange inter- 
action. A large variety of GW implementations exist, in 
which quite different levels of approximation are used. 
We will explain our procedure here and benchmark it on 
simple hydrides in the next section. 

The GqWq approximation is defined by constructing 
the self-energy £ from the orbitals and eigenvalues ob- 
tained in a DFT calculation within the local density ap- 
proximation (LDA). If the QP and LDA wave functions 
do not differ significantly, i.e., ipnk( T ) ^ Vnk A ( r )> tnen 



Eq. JIJ is approximated b} 
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where h nk = (iPnk A \ 

S„ k = (V£k A l s IV^k A >- E q- C9 is non-linear in e„ k and 
it is solved by a root-searching technique^ GW calcu- 
lations are computationally demanding, so pseudopoten- 
tials are used to represent the ion cores and only the va- 
lence electrons are treated explicitly. Calculations within 
this scheme have been applied to a wide range of semicon- 
ductors and insulators*^ They lead to band gaps that are 
usually within 10% of the experimental values, although 
occasionally somewhat larger deviations are founds 

In principle, an overlap between core and valence 
charge densities contributes to the screening, and thus 
to the self-energy. This contribution is neglected in a 
pseudopotential approach, but a simple estimate of its 
effect is made by adding 



(V xc [Pv + Pc] - V xc [p v ]) r 



(3) 



to the QP energies, where V xc is the LDA exchange- 
correlation potential, p c<v are the core and valence charge 
densities, and nk indicates the expectation value with re- 
spect to an LDA wave function as in Eq. @42i2i 

The QP equation, Eq. (QJ, is not related to DFT, but 
the scheme outlined above depends on LDA eigenval- 
ues and wave functions through the GqWo approxima- 
tion for the self-energy and through the approximation 
represented by Eq. (J2J. The dependence on LDA eigen- 
values can be avoided by constructing the self-energy 
(GW) from QP energies and solve the QP equation 
self-consistently. We have previously observed that self- 
consistency on the eigenvalues is in fact vital to obtain 
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good results for small band gap semiconductors that are 
incorrectly described by LDA as being metallic42*^i For 
large band gap materials, however, this self-consistency 
does not improve upon the Go Wo results^ The depen- 
dence on LDA wave functions can be relaxed by solving 
Eq. (JJJ instead of Eq. J5J. However, for large band gap 
materials this changes the results only marginally. 42 The 
dependence of the self-energy on the LDA wave func- 
tions is not that easily avoided. Self-consistency applied 
to QP wave functions worsens the results as compared to 
the GqWo approximation^ 

Since all hydrides considered in this paper turn out 
to be large band gap insulators, we use the Go Wo ap- 
proximation to calculate the QP spectrum. Starting 
from the optimized structures we generate wave func- 
tions and eigenvalues from an LDA calculation with 
norm-conserving pseudopotentials and a plane wave ki- 
netic energy cutoff of 748 eV4£ The self-energy Go Wo 
is calculated using the real space, imaginary time 
formarismi 41 i 44 i 45 i 46 We include 350 LDA states, use a 
real space grid mesh with a typical spacing of 0.3-0.4 
A and an interaction cell parameter of ~ 25 A . The 
QP equation is solved in the approximation represented 
by Eq. We estimate that with these parameters QP 
band gaps are converged numerically to within ±0.02 eV. 



C. Macroscopic dielectric function 

Optical excitations are two-particle excitations, but ne- 
glecting excitonic effects they can be approximated by 
transitions between single particle states. There is no ex- 
perimental indication of strong excitonic effects in metal 
hydrides. For instance, in LiH the binding energy of the 
lowest lying exciton is less than 0.05 eVi^L 4 ^ In this paper 
we assume that such excitonic effects can be neglected. 
In addition we neglect local field effects. 

If we consider quasi-particles as independent particles, 
then the imaginary part of the macroscopic dielectric 
function obtains the simple forr" ' 
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where q gives a direction, wk (ck) label single parti- 
cle states that are occupied (unoccupied) in the ground 
state, and e, u are the single particle energies and the 
translationally invariant parts of the wave functions; V 
is the volume of the unit cell. We have assumed spin 
degeneracy. 

Almost all optical data on hydrides are obtained on 
micro- or nano-crystalline samples whose crystallites 
have a significant spread in orientation. The most rele- 
vant quantity then is the directionally averaged dielectric 
function, i.e., e^ 2 \uj) averaged over q. Eq. involves 
calculating u c k+ q for small q and each k and extrapolat- 
ing to q = 0. Details can be found in Ref. |5^ 



The summation over the Brillouin zone in Eq. PJ is 
performed using a weighted tetrahedron schemed We 
found that this scheme allows for a faster convergence 
with respect to the number of k-points than various 
smearing methods. To calculate the dielectric tensor, we 
use the same plane wave kinetic energy cutoff and k-point 
mesh as for the DFT/GGA calculations. The number of 
empty bands included is sufficiently large as to describe 
all transitions up to at least 16 eV. 

If the imaginary part of the dielectric function, (uj), 
is calculated for all frequencies u>, then the real part, 
e:W can be obtained by a Kramers-Kronig transform. 
The static component eW(0) = e°° can also be calculated 
using density functional perturbation theory. Since the 
latter calculation includes local field effects, comparing 
e°° obtained with the two techniques is a way of assessing 
the importance of these effects^ 

In order to produce the right band gap one should use 
the GW QP energies in Eq. QJ. The k-point mesh used 
in an ordinary GW calculation is not sufficiently dense 
to obtain an accurate dielectric function, however, and it 
is computationally very expensive to increase the density 
of this mesh. As we will show below, the main difference 
between the GW and the GGA energies for the systems 
studied, is the size of the band gap between the occupied 
and the unoccupied states. The differences between the 
dispersions of the GW and GGA bands are relatively 
small. Therefore we adopt a simple "scissors" operator 
approximation for the energy differences in Eq. Q, 51 
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III. RESULTS 
A. Simple hydrides 

1. Test calculations 

Calculations on LiH and NaH are relatively straight- 
forward because these compounds have a simple rocksalt 
structure. They can be used to benchmark the calcula- 
tions. The calculated single particle band gaps of LiH 
and NaH are listed in Table HJ As is usual, DFT severely 
underestimates the gap, with LDA giving smaller values 
than GGA. Our calculated GW gaps for LiH and NaH are 
close to those obtained in recent PAW all-electron GW 
calculationsi 54 ' 55 As stated in the previous section, our 
calculations use pseudopotentials and take into account 
explicitly the valence electrons only. The GW band gaps 
as calculated from the QP energies obtained by solving 
Eq. @, are indicated by Ef w in Table [Q They are 
somewhat larger than the PAW values. If one applies 
the core correction of Eq. our GW band gaps be- 
come somewhat smaller than the PAW values, as shown 
by the column marked Ef w > cmc in Table |T| The differ- 
ences between the pseudopotential and the PAW GW 
gaps are small, however, i.e. of the order of 2-4% and 
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TABLE I: Single particle band gaps E g (eV) of simple hy- 
drides from DFT (GGA and LDA) and GW calculations. 
£,Gw,corc re f ers to a ppiyi n g the correction of Eq. (J^J . 





r^GGA 


ttiLDA 
h 9 






GW lit. 


LiH 


3.00 


2.61 


4.75 


4.54 


4.64", 5.24", 5.37 c 












4.99 oxp,d 


NaH 


3.79 


3.42 


5.87 


5.50 


5.68 c 


MgH 2 


3.79 


3.36 


5.64 


5.32 


5.25 ; , 5.58 9 


A1H 3 


2.18 


1.79 


4.31 


3.54 





°Ref. |H 

6 COHSEX, Ref. EH 
c Ref. 

^Experimental gap at T = 4.2 K, Ref. 

e Ref. |H 

■^cited in Ref. [57] 

»Ref. |H 



the PAW values are in between the E^ W ' COIC and Ef w 
values. 

The band gap of LiH given in Ref. |48| has been calcu- 
lated using the rather crude COSHEX approximation, 
which is known to lead to a much higher gap^ The 
value given in Ref. |56j is much higher than that obtained 
in other GW calculations, including ours, for which we 
have no explanation. The experimental band gap of LiH 
is higher than the calculated GW values^ but the dif- 
ference is within the usual 5-10%i^2i^ To our knowledge 
no experimental data are available for NaH. 

The band gap of MgH 2 is calculated for the optimized 
rutile or a-structure^ which is the most stable struc- 
ture at room temperature and ambient pressured The 
£GW,core vame we obtain is very close to that cited in 
Ref. A recently obtained PAW value for the band 
gap of MgH 2 is again in between our E^ W ' COTC and Ef w 
values. 58 For AIH3 no other data are available to our 
knowledge. One can observe that core effects are rela- 
tively large in this compound. 

In the following we will use the Eg W,core values. The 
validity of using the scissors operator approximation for 
calculating the optical response, see Eq. [5] is illustrated 
by comparing band widths calculated with DFT and 
GW. The valence band widths are given in Table [H] 
The difference between the GW and the LDA values is 
within 3% and the difference between the GW and the 
GGA values is within 10%. Note that the latter is on 
the same scale as the difference between the LDA and 
the GGA values. This accuracy is acceptable for our 
purposes. We have also checked in more detail that the 
dispersions of the individual bands in the DFT and GW 
band structures are very similar. 

Table llD also lists the static dielectric constant calcu- 
lated with and without local field effects. The small dif- 
ferences between these numbers indicate that it is rea- 
sonable to neglect local field effects in calculating the 
dielectric function. 



TABLE II: Valence band widths (eV) from DFT (GGA and 
LDA) and GW calculations. The directionally averaged real 
part of the static dielectric constant, calculated with (LFE) 
and without local field effects. 





GGA 


LDA 


GW 


00 

e LFE 


e°° 


LiH 


5.41 


5.62 


5.81 


4.28 


4.34 


NaH 


3.58 


4.00 


3.99 


3.03 


3.06 


MgH 2 


6.34 


6.62 


6.66 


3.90 


3.98 


A1H 3 


8.60 


8.82 


8.92 


4.43 


4.55 



2. LiH and NaH 

The calculated optimized lattice parameters of LiH and 
NaH in the rock salt structure are 4.02 and 4.83 A, re- 
spectively. These values are somewhat smaller than the 
experimental lattice parameters of 4.09 and 4.91 A due to 
neglecting the zero point motions of the hydrogen atoms, 
as discussed in Ref. 

The band structures and the directionally averaged 
£^ 2 \uj) of LiH and NaH are shown in Fig. ^ The va- 
lence bands in both LiH and NaH are strongly dominated 
by hydrogen, which reflects the ionic character of the 
bondingi^S* 5 ^ 5 ^* 5 ^ The conduction bands have a mixed 
hydrogen and cation character. The lattice parameter of 
LiH is significantly smaller than that of NaH. Because 
of the smaller distance between the hydrogen atoms the 
band widths in LiH are generally larger than in NaH. 
This is most easily observed in the valence band, whose 
dispersion is quite similar in LiH and NaH, but the {GW) 
valence band width in LiH is 5.81 eV, whereas in NaH it 
is 3.99 eV. 

The conduction bands of the two compounds are qual- 
itatively different. In LiH the conduction band minimum 
is at X, whereas in NaH it is at L, which causes the gap 
in LiH to be direct, whereas in NaH it is indirect. In NaH 
there is little participation of cation states in the lower 
lying conduction bands. A calculation on an fee lattice 
of H~ ions with the NaH lattice parameter and a homo- 
geneous background charge instead of the Na + cations, 
gives essentially the same band structure. In LiH the 
participation of the Li cations to the conduction band is 
larger. The conduction band minimum at X is lowered 
in energy because here the Li-2p state participates in a 
bonding combination with hydrogen states. Similarly, 
the conduction band minimum at L is raised in energy 
because the Li-2s state contributes to an anti-bonding 
combination with hydrogen states. 

The differences between the band structures of LiH 
and NaH lead to markedly different dielectric functions, 
as is shown in the lower panel of Fig. ^ The calculated 
direct optical gap in LiH is 4.54 eV, whereas in NaH it 
is a much larger 6.37 eV. The main peak in the NaH 
spectrum originates from transitions at L and X, whose 
onset is at comparable photon energies. Both valence and 
conduction bands of NaH have hydrogen character and 
the oscillator strength of these transitions is rather large. 
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FIG. 1: Upper panels: electronic band structures of LiH and 
NaH. The zero of the energy scale is at the top of the valence 
band. Lower panels: imaginary parts of the directionally av- 
eraged macroscopic dielectric functions of LiH and NaH. The 
calculated optical gaps of LiH and NaH are 4.54 and 6.37 eV, 
respectively. Unless explicitly stated otherwise, the results 
presented in the figures are based upon GGA calculations 
modified by a scissors operator extracted from the GW re- 
sults. 
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FIG. 2: Upper panels: electronic band structures of MgH 2 
and AIH3. The zero of the energy scale is at the top of the 
valence band. Pi and P2 correspond to the points (0.5,-1,0.5) 
and (0.5,0.5,0.5). Lower panels: imaginary parts of the direc- 
tionally averaged macroscopic dielectric functions of MgH 2 
and AIH3. The calculated optical gaps of MgH 2 and AIH3 
are 6.19 and 3.54 eV, respectively. 



The result is a sharp and dominant peak just above the 
onset of the optically allowed transitions. 

The optical spectrum of LiH has more structure. The 
onset at 4.54 eV is due to transitions at X, and at a 
somewhat higher energy transitions near K and W con- 
tribute. At energies ~ 9 eV transitions at L become 
allowed, which results in a peak in the dielectric function 
at that energy. Since the bands of LiH have a larger dis- 
persion than those of NaH, the optical response of LiH 
is spread out over a larger energy range. 



3. MgH 2 

a-MgH 2 has the rutile structure, i.e. space group 
P42/mnm with Mg, H atoms in 2a, 4/ Wyckoff positions, 
respectively, and two formula units per unit cell. The op- 
timized calculated lattice parameters are a — 4.52 A and 
c = 3.01 A, with the H atoms at x = 0.304. This is in 
good agreement with the experimental values a = 4.50 A, 
c = 3.01 A and x = 0.304i£i The magnesium atoms are 
sixfold octahedrally coordinated by the hydrogen atoms 
at distances between 1.94 and 1.96 A. Each MgH 6 oc- 
tahedron shares the hydrogen atoms at its corners with 
neighboring octahedra. Each hydrogen atom is shared 
by three octahedra and is therefore coordinated by three 
magnesium atoms. 



The band structure of MgH 2 is shown in Fig. [21 Our 
results are in general agreement with those obtained in 
previous workiSiS There is a small hybridization be- 
tween the H and the Mg states in the valence bands, 
with the lowest two valence bands having some Mg s and 
the highest two having some Mg p character, respectively. 
As in the case of LiH and NaH, the valence bands have 
a dominant hydrogen character, however. The conduc- 
tion bands have a mixed magnesium hydrogen character 
and the bottom of the conduction band has a consider- 
able Mg 3s contribution. MgH 2 has a calculated indirect 
gap of 5.32 eV, see Table whereas the direct optical 
gap is 6.19 eV. The experimental optical gap obtained 
in Ref. |S3 is 5.6 ± 0.1 eV, which would indicate that our 
GW result overestimates the gap by 10%. 

The lower panel of Fig. [2] shows the dielectric func- 
tion of MgH 2 . The onset of optical transitions occurs 
almost at the same energy at various regions through- 
out the Brillouin zone, which results in a steep rise of 
the dielectric function and a peak close to 8 eV. At ener- 
gies above 9 eV transitions from the lower valence bands 
start to play a role, see for instance the interval Z-M in 
the upper panel of Fig. [21 This results in a shoulder in 
the dielectric function at ~ 9 eV. Finally, the shoulder at 
~ 11 eV in the spectrum involves transitions to higher ly- 
ing conduction bands, associated with rather delocalized 
states having considerable Mg character. 
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I AIH 3 

AIH3 has a rhombohedral structure with space group 
R3c and the Al, H atoms in the 6b, 18e Wyckoff positions, 
respectively. The unit cell contains two formula units. 
The optimized lattice parameters are a — 4.49 A and 
c = 11.82 A, with the H atoms at x = 0.623. These values 
are in good agreement with the experimental values a — 
4.45 A, c = 11.80 A and x = 0.628^ The interatomic 
Al-Al distances in the ab plane and along the c axis are 
4.45 A and 3.24 A, respectively. The aluminium atoms 
form a distorted face centered structure, where each Al 
atom is octahedrally coordinated by H atoms with an Al- 
H distance of 1.75 A. Each AlB^ octahedron shares its 
corners with neighboring octahcdra and each H atom at a 
corner forms a bridge between two Al atoms. Since these 
bridges are not linear, i.e. the Al-H-Al bond angle is 141°, 
the octahedra are tilted with respect to one another. 

The band structure of AIH3 is shown in Fig. There 
is hybridization between H and Al states, but the six 
valence bands are dominated by hydrogen, as are the 
valence bands of the other hydrides. In contrast to MgH 2 , 
AIH3 has a direct band gap, which is located at T. The 
band gap is 3.54 eV, which is notably smaller than the 
gap in the other hydrides discussed above. This is caused 
by a single conduction band that disperses to ~ 2 eV 
below the other conduction bands. This band has a large 
Al 3s contribution. 

The dielectric function of AIH3 is shown in the lower 
panel of Fig. [21 Although the optical response starts at 
the direct gap of 3.54 eV, see the inset in Fig. [21 it reaches 
significant values only above 6 eV. The weak response 
between 3.54 and 6 eV is caused by the fact that only a 
single conduction band contributes with a low density of 
states. Moreover, since that band has Al 3s character, 
whereas the valence bands have dominant H character, 
the oscillator strength of these transitions is small. The 
dielectric function rises sharply above 6 eV and peaks 
above 7 eV. The spectrum has a distinct broad shoulder 
between 9 and 10 eV and also some weaker shoulders at 
higher energies. 

In order to interpret the dielectric function it is instruc- 
tive to analyze the total density of states (DOS) and the 
local density of states (LDOS) of AIH3. These are shown 
in Fig. The valence DOS has a sharp peak just below 
the Fermi level and two broad peaks at ~ —2.5 eV and 
~ —6 eV below the Fermi level. These three peaks orig- 
inate from respectively 3s, 3p and 3d aluminium states 
hybridizing with the Is hydrogen states, as can be ob- 
served in the upper two panels of Fig. [21 In the solid the 
states are broadened into strongly overlapping bands. In 
the dielectric function transitions from the highest two 
of these valence peaks gives rise to the structure between 
6 and 10 eV in Fig. [51 The energy associated with tran- 
sitions from the lowest valence peak is too high to give 
any significant contribution to the dielectric function. 

Qualitatively the optical spectra of AIH3 and MgH 2 
show some similarity, despite the difference in structure 



A1H 3 




E - Ep (eV) 

FIG. 3: (Color online) Local densities of states (LDOS) in 
atomic angular momentum projection; bottom panel: total 
density of states. The Fermi energy is at the top of the valence 
band. 

between these materials. In both cases the dielectric 
function sharply rises above 6 eV and peaks between 7 
and 8 eV. The spectrum of AIH3 is broader due to a larger 
dispersion of the bands, reflecting the somewhat denser 
packing of the hydrogen atoms in this compound. From 
Fig. ^ one observes that the dielectric function of NaH 
also rises steeply above 6 eV and peaks just above 7 eV. 
The spectrum of NaH is narrower than that of MgH 2 and 
AIH3 , reflecting the less dense packing of hydrogen atoms 
in this compound, which results in a smaller band disper- 
sion. Only the spectrum of LiH is qualitatively different 
as it rises below 5 eV in a broad shoulder. As discussed 
in Sec. lIII A~2l there is a significant contribution from the 
Li states in this case. 

B. The binary hydrides LigAlHe and Na 3 AlH6 

The optimized atomic positions of LisAlHg and 
Na 3 AlH 6 are given in Table [nil Both Li 3 AlH 6 and 
Na3AlH6 consist of a stacking of AlHg octahedra and 
alkali cations. The octahedra are slightly distorted with 
Al-H distances of 1.75 A in Li 3 AlH 6 and 1.78 to 1.80 A 
in NasAlHg. These compounds contain a relatively large 
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TABLE III: Optimized atomic positions in the binary hy- 
drides. The labels "la" etc. refer to the Wyckoff positions. 
The cell parameters are taken from the references. The struc- 
tures are in good agreement with previous experimental and 
theoretical work, 8 i 26 i 27 . 28 . 64 i 65 i 66 i 67 . 68 . 69 



TABLE IV: Single particle band gaps E g (eV) of the binary 
hydrides calculated with DFT (GGA and LDA) and GW. 
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fraction of alkali cations. Since sodium atoms are larger 
than lithium atoms, the distance between the AlHg oc- 
tahedra in NasAlHg is considerably larger. The Al- 
Al distance in NaaAlHg is 5.59 A, whereas in LisAlHg 
it is 4.88 A. As for the simple hydrides discussed in 
Sec. IIII A~2l this size effect leads to noticeable differences 
in the electronic structure and the optical properties of 
Li 3 AlH 6 and Na 3 AlH 6 . 

The local electronic densities of states (LDOS) of both 
compounds is given in Fig. As in AIH3 the valence 
bands have dominant hydrogen character, although there 
is Al character mixed in. The splitting into three peaks 
with approximate relative intensity 1:3:2 can be inter- 
preted in terms of an octahedral ligand field splitting. 28 
The peaks correspond respectively to the 3s, 3p and 
3d(e g ) aluminum states hybridizing with the Is hydro- 
gen levels in the AlHg octahedra. As can be observed 
in Fig. 0] the spd splitting in LisAlHg and NasAlHg is 
comparable, which reflects the similarity of the AlHg oc- 
tahedral structure in both compounds. Comparing to 
Fig. |3 one observes that the splitting is also comparable 
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to that in AIH3, again suggesting the similarity in the 
octahedral structure. 

The interaction between the octahedra in the solid 
results in a broadening of the three peaks. Unlike in 
AIH3 the A1H6 octahedra in LiaAlHg and NaaAlHg are 
not directly connected, which limits the interaction and 
the broadening. Therefore, the three peaks remain non- 
overlapping. One expects their widths to increase as the 
distance between the octahedra decreases and indeed the 
valence peaks in the DOS of Li 3 AlH 6 are wider than in 
NasAIHg. The conduction band of both compounds is 
rather featureless up to at least 22 eV. There is somewhat 
more cation s character at the bottom of the conduction 
band in NasAlHg, whereas the bonding in LiaAlHg prob- 
ably has a somewhat more covalent character, as in the 
simple hydrides. 

The calculated band gaps are given in Table ITVI As 
discussed in Sec. IIII A II DFT calculations severely un- 
derestimate the gap, with LDA giving a ~ 0.5 eV smaller- 
value than GGA. The most important results are in the 
last two columns of Table II VI which give the GW single 
particle gap and the direct optical gap. The single par- 
ticle gap is indirect in L13AIH6 and direct in Na3AlH6. 73 
The fairly large difference between the gaps of Li 3 AlH 6 
and NaaAlHg is striking. Moreover, the fact that the 
gap of LiaAlHg is larger is somewhat counterintuitive. 
Naively one would expect that the larger broadening of 
the A1H 6 octahedron levels discussed above would nar- 
row the gap, since it leads to a larger valence band width. 
The origin of the band gap difference between L13AIH6 
and Na3 A1H 6 is discussed in Sec. IIVI 

The directionally averaged dielectric functions of 
LisAlHg and NasAlHg are shown in Fig. [5] The dielectric 
function of L13AIH6 has a shoulder starting above 6 eV, 
a peak just above 8 eV and a shoulder below 12 eV. Since 
the conduction band DOS is rather uniform and feature- 
less up to at least 22 eV, these features in the dielectric 
function can be directly linked to transitions from the 
three octahedron valence peaks. 

Despite the much smaller band gap of NasAlHg the di- 
electric function starts to increase appreciably only at an 
energy between 5 and 6 eV, which is not that much lower 
than in LisAlHg . Transitions from the top two valence 
peaks give rise to the complicated pattern between 6 and 
9 eV; transitions from the third valence peak gives the 
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FIG. 4: (Color online) Local densities of states (LDOS) of LisAIHe and Na3 Allied- For clarity the area under the d-line is 
shaded. The Fermi energies are at the top of the valence band. The conduction band DOS is almost constant up to at least 
22 eV. 
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FIG. 5: Imaginary parts of the directionally averaged macro- 
scopic dielectric functions of LiaAIHe and Na3AlH 6 . 



above 10 eV. Qualitatively these spectra have a resem- 
blance to that of AIH3, see Fig- El reflecting the dominant 
role played by the AlHg octahedra. 



C. The binary hydrides LiAlH 4 , NaAlH 4 and 
Mg(AlH 4 ) 2 

The optimized structures of LiAlH4, NaAlfi4 and 
Mg(AlH 4 ) 2 are given in Table Hill For LiAlH 4 and 
NaAlH4 we have used the experimental unit cells and 
optimized the atomic positions only; for Mg(AlH 4 ) 2 we 
have also optimized the size and shape of the unit celli 
All three materials consist of a packing of AIH4 tetrahe- 
dra and alkali or alkaline earth cations. The tetrahedra 
are slightly distorted and the Al-H distances vary from 
1.62 to 1.65 A in LiAlH 4 , 1.64 A in NaAlH 4 and from 1.60 
to 1.62 A in Mg(AlH 4 ) 2 . Unlike the two compounds dis- 
cussed in the previous section, the volume fraction taken 
up by the cations is relatively small and the distance be- 
tween the AIH4 tetrahedra is hardly influenced by the size 
of the cations. The Al-Al distance is 3.75 A in L1AIH4, 
3.78 A in NaAlH 4 and 3.86 A in Mg(AlH 4 ) 2 . 

The LDOS of LiAlH 4 , NaAlH 4 and Mg(AlH 4 ) 2 is 
shown in Fig. As for the compounds discussed be- 
fore, the valence bands have dominant hydrogen char- 
acter with some Al character mixed in. The splitting 
into two peaks with approximate relative intensity 1:3 is 
due to a tetrahedral ligand field splitting of the Al 3p 
and 3s levels hybridized with H Is levels in AIH4. 28 The 
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FIG. 6: (Color online) Local densities of states (LDOS) of LiAlH 4 , NaAlH 4 and Mg(AlH 4 ) 2 !£- For clarity the area under the 
d-line is shaded. The Fermi energies are at the top op of the valence band. The conduction band DOS is almost constant up 
to at least 22 eV. 



splitting is comparable in all three compounds, which re- 
flects the similarity of the structure and bonding of the 
AIH4 tctrahedra in these compounds. The interaction 
between the tetrahedra causes a broadening of these lev- 
els in the solid. Fig.[S]shows that also the resulting band 
widths of these valence states is comparable in all three 
compounds. Apparently the widths are not extremely 
sensitive to the details of the structure, which are quite 
different for LiAlH 4 , NaAlH 4 and Mg(AlH 4 ) 2 . They are 
sensitive to the distance between the tctrahedra, but this 
is comparable for the three compounds. 

Compared to the valence bands, the features in the 
conduction bands are less distinct. Both in NaAlH 4 and 
in Mg(AlH 4 ) 2 the bottom of the conduction band has 
considerable s character derived from the empty 3s state 
of the cation. The conduction band in LiAlH 4 is fea- 
tureless. The LDOS on the Al atoms is very similar in 
NaAlH 4 and in Mg(AlH 4 ) 2 , but there are small differ- 
ences with LiAlH 4 . There is significant Al d character 
in the valence band in the latter compound, and almost 
none in the other two compounds. In the conduction 
band of LiAlH 4 there is a considerable Al s contribution, 
and much less in the other compounds. 

In order to evaluate these differences we have also cal- 
culated the LDOS for the LiAlH 4 structure with the 
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FIG. 7: Imaginary parts of the directionally averaged 
macroscopic dielectric functions of LiAlH 4 , NaaAIHe and 
Mg(AlHJ 2 . 

Li + ions replaced by a uniform positive background. The 
LDOS on the Al atoms then becomes very similar to that 
in the NaAlH 4 and Mg(AlH 4 ) 2 compounds. All these fea- 
tures indicate that NaAlH 4 and Mg(AlH 4 ) 2 can be con- 
sidered as ionic compounds, i.e., as a packing of A1HJ an- 
ions and Na + or Mg 2+ cations, whereas in LiAlH 4 there 
may be a stronger covalent contribution. 

The calculated band gaps are given in Table ITVl The 



10 



GW single particle band gaps of the three compounds 
are almost the same and also in the optical gaps there 
is very little differenced This similarity indicates that 
the electronic structure around the band gap is foremost 
determined by the (AIH4) tetrahedra. The distances 
between these tetrahedra are similar in these three com- 
pounds and apparently the detailed differences in their 
packing are relatively unimportant. 

This conclusion is strengthened by the dielectric func- 
tion, which is shown in Fig. The maximal dielectric 
response of Mg(AlH 4 ) 2 is somewhat smaller than that of 
LiAlH4 and NaAlH4, but the shape of the three curves 
is remarkably similar. The double peak structure of the 
valence band of the LDOS, which appears in all three 
compounds in Fig. El is almost washed out in the dielec- 
tric response. Transitions from the lowest valence band 
can be recognized only as a faint shoulder near 10 eV. 
It is then not surprising that the smaller differences be- 
tween the LDOS of LiAlEU and the other two compounds 
do not influence the dielectric functions much. 



IV. DISCUSSION 

The electronic structure and the dielectric function of 
the binary compounds discussed in Sees. IIII 51 and ITTTCl 
are foremost determined by the lattice of (AlHg) - octa- 
hedra and (AIH4) tetrahedra, respectively, whereas the 
cations have a minor influence. In this section we will dis- 
cuss this proposition in more detail. The exact value of 
the band gap is not important in this discussion, only its 
relative variation with structure and composition. Since 
the latter is described qualitatively by DFT/GGA cal- 
culations, see Table llVl we will only use GGA results in 
this section. 

We have calculated the dielectric functions of lat- 
tices of (A1H 4 )~ tetrahedra in the LiAlH 4 , NaAlH 4 
and Mg(AlH 4 ) 2 structures, but with the Li + , Na + and 
Mg 2+ cations replaced by a uniform positive background 
charge. The results are compared to the dielectric func- 
tions of the real compounds in Fig. |SJ It can be observed 
that removing the cations in LiAlELi hardly changes the 
dielectric function. Removing the cations in NaAlH4 re- 
sults in a slight shift of the dielectric response to higher 
energies. This is related to the disappearance of the 
peak at the bottom of the conduction band, which has 
a sodium s character, see the middle panel of Fig. 
Also in Mg(AlH 4 ) 2 removing the cations results in small 
changes in the dielectric function only. These are mainly 
caused by the disappearance of the magnesium related 
peaks at the bottom of the conduction band and a re- 
sulting flattening of the conduction bands, see the right 
panel of Fig. El 

In conclusion, although removing the cations results 
in small changes in the conduction band, overall the di- 
electric function changes very little, which means that it 
is foremost determined by the lattice of (AIH4) anions. 
Such a behavior is not uncommon for ionic compounds. 
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FIG. 8: (Color online) Dielectric functions of MAIH4 com- 
pounds (black, solid) and of the corresponding systems with 
Li + , Na + , and Mg 2 + ions substituted by a uniform back- 
ground charge (red, dashed). These results are based upon 
GGA calculations without scissors operator cerrection. 



For instance, in alkali halides such as NaCl both the top 
of the valence band and the bottom of the conduction 
band are determined by the anion latticeiS^ 

Going one step further one can correlate the relative 
size of the band gap with the distance between the an- 
ions. We will illustrate this using the DOS of Mg(AlH 4 ) 2 , 
which is shown in Fig.^Ja). Replacing the Mg 2+ ions by 
a uniform positive background does not change the DOS 
significantly, as can be observed in Fig.^b). In Fig.EJc)- 
(e) the cell parameters are increased while the geometry 
of the (AIH4) tetrahedra is fixed. As the distance be- 
tween the anions increases, the band widths of all bands 
decreases, but those of the valence bands decrease much 
more rapidly. At a large distance the DOS is essentially 
that of an isolated (A1H4) - tetrahedron, where the va- 
lence states are split into states of s and p symmetry due 
to the tetrahedral ligand field. The upper states of p sym- 
metry show a small splitting, since the tetrahedron has 
a small trigonal distortion. The gap between the high- 
est valence state and the lowest conduction state in the 
isolated tetrahedron is only ~ 0.5 times the band gap in 
the solid, compare Figs.OJb) and^e). 

In general, the larger the distance between the anions, 
the smaller the gap. This result seems to be somewhat 
counterintuitive, as at the same time the valence and con- 
duction band widths are decreasing and in general this 
would increase the gap. The result can be understood 
in terms of the electrostatic (Madelung) potential^ De- 
creasing the lattice constant, the electrostatic potential 
on the anions becomes more attractive to electrons, due 
to the closer packing of the cations. The same argu- 
ment also holds for a uniform positive background, i.e., 
a decreasing lattice constant leads to a more attractive 
Madelung potential on the anions. All states on the an- 
ions experience this potential and lower their energy. The 
size of this shift, however, depends upon the degree of lo- 
calization of the state. If a states is completely localized, 
the Madelung potential will maximally lower its energy. 
On the other hand, if a state is completely delocalized, 
its energy shift is zero, since the system as a whole is 
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FIG. 9: GGA densities of state of (a) Mg(AlH 4 ) 2 and (b) 
with the Mg 2+ ions replaced by a uniform background charge. 
While fixing the geometry of the (AlHU) - ions, the lattice is 
expanded by a factor of 1,5 (c), 2 (d), and 5 (e). 



charge neutral. 

The key point is that in our systems the valence states 
are much more localized than the conduction states. This 
is immediately evident from Fig. GHb)-(e), where the va- 
lence band widths decrease much faster with an increas- 
ing lattice constant than the conduction band widths. 
As a result, the more localized valence states increase 
their energy significantly faster with an increasing lattice 
constant than the conduction band states. Since this ef- 
fect is much larger than the effect of the decreasing band 
widths, increasing the lattice constant results in a smaller 
band gap. 

This is quantified in Fig. ^3 where the band gaps of 
the 2, 3, 4, 5 and 6 times expanded lattice are fitted to 
an expression Ae 2 /d, with d the distance between the 
anions. The constant A = 1.57 represents the difference 
in localization of the valence and conduction states. This 
simple model breaks down if the localization of the states 
strongly depends upon d, i.e. in the 1.5 and 2 times 
expanded lattice. It occurs if d becomes sufficiently small, 
see Fig.EDV). 

This concept can be used to interpret the trend in 
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FIG. 10: Calculated GGA band gaps of the expanded 
(AIH4) "-lattice (circles) fitted to Ae 2 jd (dashed line), whith 
d the distance between the anions; A — 1.57. 



the band gaps of the binary compounds, see Table ITVl 
In Li AIH4 , NaAlH 4 and Mg(AlH 4 ) 2 the Al-Al distance, 
which is a measure for the distance between the anions, 
is almost the same, so their band gaps are very close. 
The Al-Al distance in LisAlHg is much smaller than in 
NasAlHg, which explains the larger band gap in the for- 
mer compound. 



V. SUMMARY 

In this paper the electronic structures and dielectric 
functions of the simple hydrides LiH, NaH, MgH 2 and 
AIH3, and the complex hydrides LisAlHg, NaaAlHg, 
LiAlH 4 , NaAlH 4 and Mg(AlH 4 ) 2 , have been studied by 
first principles calculations. The equilibrium structures 
of these compounds are obtained from DFT/GGA total 
energy minimizations. GW calculations within the QP 
approximation provide the single particle excitation en- 
ergies, i.e., the electronic band structures. We use the 
GqWq approximation based upon LDA wave functions 
and eigenvalues. The difference between the dispersions 
of the GW and the GGA bands is less than 10%. There- 
fore, the band structures are well represented by GGA 
band structures that are corrected by applying a scissors 
operator between occupied and unoccupied states in or- 
der to obtain the GW band gap. From the single particle 
wave functions we then calculate the directionally aver- 
aged dielectric functions within RPA, neglecting exciton 
effects. We also neglect local field effects, but from cal- 
culations on static dielectric constants we conclude that 
this is a reasonable approximation. 

All compounds are large gap insulators with band gaps 
that vary from 3.5 eV in A1H 3 to 6.5 eV in the MAIH4 
compounds. In all cases the valence bands are domi- 
nated by the hydrogen atoms, whereas the conduction 
bands have mixed contributions from hydrogen and metal 
cation states. The band gap in LiH, AIH3 and NasAlHg 
is direct, whereas in all the other compounds it is indi- 
rect. The optical response of most compounds is qual- 
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itatively similar, notwithstanding sizeable differences in 
their band structure and band gap. The dielectric func- 
tion (w) rises sharply at photon energies correspond- 
ing to ~ 6 eV, and around ~ 8 eV it has a strong peak 
reaching values in the range 10-15. In the direct gap ma- 
terials £^{lo) has a weak tail going to lower energies. 
Between <~ 8 and ~ 12 eV, e^(uj) gradually decreases 
to a value < 2 at 12 eV. Most of the materials specific 
optical information can be found in this energy range, 
albeit in the form of relatively weak shoulders in e^(uj). 

The electronic structure and the optical properties of 
the aluminium compounds can be interpreted in terms 
of aluminium hydride complexes, i.e., AlHg octahcdra 
in AIH3, LisAlHg and NasAlHg, and AIH4 tetrahcdra 
in L1AIH4, NaAlH4. Explicit calculations on lattices of 
these complexes, without the Li, Na, and Mg cations, 
show that the latter have a relatively small effect on 
the DOS and on the optical response. The distance be- 
tween the (AIH4) - tetrahedra in L1AIH4, NaAlELt and 
Mg(AlH 4 ) 2 is almost the same. Since the interaction be- 
tween the tetrahedra is then similar, this explains why 
the optical spectra of these compounds are very similar. 

The same reasoning can be applied to Li 3 AlH 6 and 



NasAlHg in terms of a lattice of (AlHg) 3 octahedra. 
However, the distance between the octahedra is smaller 
in the Li compound because of the smaller size of the 
cation. The band gap then becomes larger, which can be 
understood from the influence of the increased Madclung 
potential on the more localized valence states. 
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